Constraining the variation of the coupling constants 
with big bang nucleosynthesis 
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We consider the possibility of the coupling constants of the 57/(3) x SU (2) x (7(1) gauge interactions 
at the time of big bang nucleosynthesis having taken different values from what we measure at 
present, and investigate the allowed difference requiring the shift in the coupling constants not 
violate the successful calculation of the primordial abundances of the light elements. We vary gauge 
couplings and Yukawa couplings (fermion masses) using a model in which their relative variations 
are governed by a single scalar field, dilaton, as found in string theory. The results include a limit 
on the fine structure constant —6.0 x 10~ 4 < Achem J '&EM < 1.5 x 10~ 4 , which is two orders stricter 
than the limit obtained by considering the variation of cxem alone. 



I. INTRODUCTION 

Big bang nucleosynthesis (BBN) is one of the most im- 
portant tools for probing the early universe. The stan- 
dard big bang nucleosynthesis (SBBN) predicts the pri- 
mordial abundances of light elements (D, 3 He, 4 He and 
7 Li) with the only free parameter, the baryon-to-photon- 
ratio r\. For r\ ~ O(10™ 10 ), they are remarkably consis- 
tent with the observations of the abundance which extend 
over 9 digits. Requiring not to vitiate this consistency, 
the stringent constraints on various theories which affect 
cosmology have been obtained. 

Among such theories constrained by BBN, there are 
some in which values of coupling constants may vary. 
The SBBN prediction assumes, along with three light 
neutrino species and no lepton asymmetry, that physical 
parameters involved in the calculation (the fine struc- 
ture constant czem, the Fermi constant Gf, the electron 
mass m ei etc.) are the same for now and the BBN time. 
However, this is not an obvious choice because there are 
numbers of ways in which they could have varied. In 
addition to the theoretical possibility, the recent analy- 
sis of quasar absorption lines found possible evidence for 
variation in a^M 0- 

In this paper, we consider a model taken from string 
theory where these coupling constants are all related to 
the expectation values of dilaton field <& and could in 
principle vary with time. For example, very general ar- 
guments indicate that the dilaton cannot be stabilized at 
a value we would characterise as corresponding to weak 
coupling Qj. Also, in certain popular models for stabilis- 
ing the dilaton using gaugino condensates, the cosmolog- 
ical evolution would almost inevitably tend to overshoot 
the desired minimum of the dilaton potential and run off 
to an anti-de Sitter vacuum ||. 

The notion that the coupling constants are determined 
by a single scalar field not only motivates the investiga- 
tion of their time variation as a probe for physics beyond 
the Standard Model, but also makes the analysis of si- 
multaneous changes in the couplings simple and concrete. 
In addition to the free parameter of the SBBN 77, we have 



only another parameter 
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tional variation of the dilaton field. After we calculate the 
primordial abundances with various values of these two 
quantities, we search a parameter region not excluded by 
the observation. Thus obtained limit on A<j>/$ is readily 
translated in the limits on the coupling constants. 

Our analysis is positioned as the extension to Ref. [Q. 
There, the dilaton dependence of the coupling constants 
is determined by the action of the heterotic string in 
the Einstein frame. The constraint is obtained from 4 He 
whose abundance and coupling dependence can be esti- 
mated without recourse to a numerical calculation. Their 
method is easy to calculate and appropriate for an order 
estimation. But we can extract more information by per- 
forming the numerical calculation, whose advantages are: 
1. the abundances of the light elements other than 4 He 
are calculated so we can use more observational data, 
especially D's, to put a constraint, 2. we can take into 
account of 77 dependence of the abundances, 3. we use re- 
alistic values for the reaction rates so we avoid the rough 
estimation of the weak reaction rate (both absolute value 
and temperature dependence) based on dimensional anal- 
ysis, 4. by Monte-Carlo simulation, we can estimate theo- 
retical uncertainty, 5. by calculating a statistical measure 
(such as x 2 ), we can objectively quantify the constraint. 
Of course, as the SBBN becomes more and more precise 
owing to the progress in the studies of nuclear reaction 
rates and of primordial abundances, numerical computa- 
tion is necessary in general to reflect the recent develop- 
ment and to produce a result close to the truth. 

In the next section, we review briefly the SBBN and 
see what kind of physical quantities are needed to make 
the prediction. In Sec. IB, we estimate their coupling de- 
pendences. In Sec. IV, we introduce a model containing 
the dilaton and we investigate the dilaton dependence of 
the coupling constants. In Sec. [V], we put a bound on 
A$/$ by calculating x 2 with observational and theoret- 
ical uncertainty from the recent data (the former very 
much dominates). In terms of the limit on Aa^M I 'ct-EMi 
we find 
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two orders more restrictive than the one found in Ref. |j| 
where only the variation of olem is taken into account. 
We discuss this originates from the unified gauge cou- 
plings, which relates electromagnetic sector to strong sec- 
tor. At the energy scale of BBN, manifest quantities are 
o-em and Aqcd which appears instead of a s t r0 ng through 
dimensional transmutation. The former linearly depend 
on the dilaton but the latter does exponentially. 4 He 
abundance is determined by the magnitude of neutron- 
proton mass difference Am which is in turn determined 
by Aqcd and this dominates the constraints by BBN. 



II. PHYSICAL QUANTITIES IN THE SBBN 

In order to see how the effects of changing coupling 
constants arise, we summarise the main points of the 
SBBN calculation. Following Ref. 0, we divide it into 
three stages. The actual calculation is performed by com- 
puter which runs the code to solve the set of ordinary 
differential equations so, of course, does not distinguish 
such stages but we see the important quantities which de- 
termine the primordial abundances of the light elements 
appear in this brief review. 

1st Stage (statistical equilibrium)^ 3> 1 MeV; (« 1 
sec) 

The energy density and the number density are dom- 
inated by the relativistic electron e~ , positron e + , neu- 
trino v, anti-neutrino v and photon 7. There are three 
types of neutrino and anti-neutrino. There is only tiny 
fraction of protons and neutrons, just ~ 10~ 9 number 
fraction. In this period, all these particles make elastic 
scattering so frequently that they are in thermal equi- 
librium and has equal temperatures. Especially, weak 
interaction cross sections are large and neutrinos are in 
thermal equilibrium with the other particles. In addi- 
tion, there are weak interaction processes interchanging 
protons and neutrons, 



n + v e «-> p + e 
n + e~"~ <-> p + v e 
n «-> p + e~ + v e , 



(2) 
(3) 
(4) 



so they are also in the chemical equilibrium. Then, the 
number ratio of proton to neutron at temperature T is 
Ha = e -{m„-m v )/T e (p, n -n P )/T w h ere rn n {m p ) is neutron 

(proton) mass and p, n (fip) is its chemical potential, but 
the last factor e^ tJ,n ~ ,J ' p ^ T is almost unity if the universe 
has no lepton asymmetry, which is the assumption of 
SBBN. So as long as equilibrium holds with tempera- 
ture T, the neutron-to-proton ratio is determined by the 
neutron-proton mass difference Am = m n — m p , 



n n _ —Am/T 



(5) 



At this stage, the nuclear reaction rates are also fast 
enough for light nuclei to be in chemical equilibrium. 
Then the abundances of the nuclear species with mass 
number A and atomic number Z is 



Y A = gAm) A - 1 n (1 - A)/2 2 {3A - 5)/2 }A 5 ' 2 ( — 

\m N 

x v A -% z Y n A - z eMB A /n 



3(A-l)/2 



(6) 



where £(3) ~ 1.202, m^ is the nucleon mass, Ba is the 
binding energy and gA is the number of degrees of free- 
dom. Since 77 is small (~ 10~ 10 ) and Ba -C T, nuclei 
abundance is negligible at this temperature. 

2nd Stage (neutron-proton freeze-out) (T ~ 0.7 MeV; 
t ~ 2 sec) 

As the universe expands and cools, the rates of inter- 
actions involving neutrinos decreases. Especially, there 
are some point when the reactions (§)-(|) practically 
stop and the chemical equilibrium break down. After 
this point, the numbers of neutrons and protons do not 
change (actually, small amount of neutrons turns into 
protons by the beta decay) and their ratio is fixed (froze 
out) at the value determined from Eq. (jSJ) with the tem- 
perature of that epoch. We call this "freeze-out temper- 
ature" and denote Tf. 

The freeze-out occurs when the rate of the reaction 
r(T) and the expansion rate of the universe H(T) be- 
comes equal: 



H(T f ) 



weak 



(Tf). 



(7) 



The expansion rate is determined by the energy density 
p through the Friedman equation, 



8irGp 



(8) 



At this stage, p is dominated by photons, neutrinos and 
relativistic electrons (T > m e = 0.511 MeV) and hence 
the cosmic density is given by p = 10.75(7r 2 /30)T 4 . The 
rate of weak interaction is 

r»eot = r(n -> p) 

= T(n + v e — > p + e~) 

+ Y(n + e + — > p + Ve) + T(n — > p + e~ + v e ). 

Using the Fermi theory of the weak interaction, this can 
be written as 



T(n p) 
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1 + e x l T » 1 + e -(x+&rn)/Te 
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We calculate the normalization factor A by using T(n — > 
p + e~ + C" e )|T=o = T n 1 where r n is the neutron lifetime. 
Integration gives A = 0.05606 r" 1 MeV -4 . From these 
expression and Eq.(^), we find Tf ~ 0.7 MeV and, from 
Eq. (|), freeze-out ratio {n n /n p ) f = e ~ Am / T / w 0.158. 

For the nuclear reaction, the situation did not change 
from the previous stage. Their abundances are still very 
small. 

3rd Stage (light- element synthesis) (0.7 MeV > T > 
0.05 MeV; 3 sec < t < 6 min) 

After the second stage, the electron-positron pair an- 
nihilation completed so there remains the photons and 
the neutrinos as relativistic particles. 

Since the early universe has lower density than inside 
of stars, there practically occurs only two body reactions. 
Thus unless the deuterons are synthesized by the reaction 
p + n — > D + 7, larger nuclei are not synthesized. This 
reaction does not proceed until the photon density is low 
enough not to photodissociate D or, in other words, the 
factor r]exp(Bi)/T) in Eq.(^) becomes O(l). It occurs at 
about T = 0.06 MeV. After that, mainly charged-particle 
reactions such asD + D^t+p and t + D — > 4 He + n 
proceed synthesizing almost all of the neutrons into 4 He 
because it has the largest binding energy among the light 
elements. Therefore, He abundance (conventionally ex- 
pressed by mass ratio) can be estimated by the frozen 
n-p ratio found at the 2nd stage: 



Yi 



He 



l + (n p /n n ) f l + e A '"/r, 



(9) 



Besides 4 He, small amount of D and 3 He and very small 
amount of 7 Li are synthesized. Their abundances are 
determined by n and the reaction rates. 

In summary, to perform BBN calculation, we need to 
input following values at BBN time: 

• neutron-proton mass difference Am 

• neutron lifetime t„ 

• nuclear reaction rates 

Mainly, the abundance of 4 He are affected by the first 
two as seen from Eqs. (0) and (|)|), and D, 3 He and 7 Li 
by the last one. 



III. THE COUPLING CONSTANTS 
DEPENDENCE OF THE BBN INPUTS 

We know the present value of the BBN input quanti- 
ties pointed out in the previous section with some exper- 
imental uncertainty. To estimate their values during the 
BBN, we try to express these quantities in terms of the 
coupling constants. 



A. neutron-proton mass difference Am 

The origin of the neutron-proton mass difference is 
traced to the electromagnetic self-energy difference and 
the d-u quark mass difference (QCD chiral symmetry 
breaking by the mass terms). The former makes proton 
heavier than neutron but the latter does the inverse and 
the total is measured accurately to be Am = m n — m p = 
1.2933318 ± 0.0000005 MeV ff. 

As argued in Ref. M the largest contribution to the 
electromagnetic part comes from the Born term of the 
Cottingham formula || and hence it can be calculated 
with relatively less uncertainty. This formula expresses 
the nucleon self-energy in the first order of czem by its 
electric and magnetic form factor. Using the formula and 
experimental data for the form factors, the difference in 
the neutron and the proton is — 0.76MeV. This value is 
proportional to (of course) olem and to 1^-qcd to have a 
proper dimension. 

On the other hand, the absolute values of u and d 
quark masses are not well known, m u = 1 ~ 5 MeV and 
ma = 3^9 MeV [|J. But knowing the electromagnetic 
contribution, the quark contribution should be 2.05 MeV. 
This is proportional to the Yukawa couplings y u ,yd and 
the Higgs expectation value (H). 

Now we can write 



Am = aa EM K QCD + b(y d - y u )(H), 



(10) 



where a and b are constants we assume not to depend 
on any coupling constants. Thus, we estimate the mass 
difference at BBN epoch Amtbu as 

Ambb„ 

= aOiEMbbnh-QCDbbn + bybbn{H)bbn 

_ _q OtEMbbn A-QCDbbn + ^ Vbbn {H)bbn J Me yj 



OtEMnow AqcD 



Unow y^O 



(ii) 



B. neutron lifetime 

Neutron /3-decay n — > p + e + 9 is very well approxi- 
mated by the one-point interaction of four particles: neu- 
tron, proton, electron and neutrino. Its coupling constant 
is denoted Gf called the Fermi coupling constant. Using 
this theory, 

t" 1 (xG 2 F J d z p e d z p v S(E e + E u + m p - m„) 

p A rn 

= G 2 F 16tt 2 / dE e E e yjEl-m 2 e (Am - E e ) 2 

= G\ \^m\ 1{ V^T(2g 4 - 9q 2 - 8) 

+ 15 g log(<7+vV-l)} (12) 
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where we defined q = ^p. Denoting the factor . .} 
/(«), 



7(g)- 1 



= (H)*(y e (H))- 5 f(q)-i = {H)-^f{ q )~\ (13) 

2 2 

The first equality follows from Gf = -jgr = ( g Jji))^ = 
jjjy? where gi is the SU(2) coupling constant and M\y is 
the weak boson mass. Therefore, we obtain 



7~n,bbn 



Ue,bbn 


-5 


(H)bbn 


-1 


' f(qbbn) ' 
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_f{Qnow)_ 



(14) 



where y e is the electron Yukawa coupling and f(q n ow) 
1.63615. 



Ref. |T(| for the latter, and they are assumed to be in- 
dependent of the couplings because there are no theoret- 
ical derivation for these. Because 3 Ke(n,p)t affects 3 He 
abundance which is not used for our analysis, our result 
does not change by neglecting its coupling dependence. 
On the other hand, 7 Be(n,p) 7 Li affects 7 Li abundance a 
lot when r\ is large so its effect may be large. However, 
since there is large observational uncertainty in 7 Li, it is 
expected that our result does not change but we have to 
know its coupling dependence to extract a reliable con- 
straint from 7 Li data when the uncertainty of 7 Li abun- 
dances decreases in the future. 



IV. THE DILATON DEPENDENCE OF THE 
COUPLING CONSTANTS 



C. charged-particle induced reaction rates 

Most of the reaction rates involved in the BBN calcu- 
lation are charged-particle induced reaction rates. Their 
clem dependence is considered in Ref. H for the rates 
discussed in Ref. (l(J. We implemented the osem depen- 
dence in the same manner updating the reaction rates re- 
cently compiled by Angulo et al. |ll[ . Since only a few re- 
action rates have resonance terms which is considered to 
be dependent of the strong coupling constant, or in other 
words since the charged-particle reaction rates influential 
for BBN are practically determined by Coulomb barrier 
penetrability, we expect that neglecting their strong cou- 
pling dependence will not affect the results a lot. 

D. neutron-induced reaction rates 

The cross section for n+p — > D+7 is calculated at ener- 
gies relevant to BBN using the effective field theory that 
describes the two-nucleon sector |l^] . The rate obtained 
by thermal averaging this theoretical cross section repro- 
duces the abundance calculation using the rate shown 
in Ref. fio|| . We exploit this theoretical formula to esti- 
mate the coupling dependence by assuming the param- 
eters which have the dimension [length]™ in the formula 
proportional to m~ n . The coupling dependence of pion 
mass m n is known from Gell-Mann-Oakes-Renner rela- 
tion to be m 2 , cx rn quar kh-QCD- The other parameters in 
the formula are the nucleon mass mjy and the deuteron 

binding energy Be- Their dependence is estimated by 

2 

m N oc Aqce and B D cx a "* d cx m quark . This B D 's 
dependence is estimated considering a simple square well 
potential for deuteron. It is also estimated by uncertainty 

principle: p ~ ~ m„, and -B D + ~ 0. 

For the other neutron-induced reactions, 3 He(n,p)t 
and 7 Be(n, p) 7 Li, we use the rate fitted from the exper- 
imental data as found in Ref. [Hi for the former and 



Next, we introduce the action which governs the varia- 
tion of the coupling constants as is used in Ref. [|J . This 
is the tree level low energy action of the heterotic string 
in the Einstein frame. 

S = J d 4 x^y^R - \d^8^ - \d^d^ 

- q,- 2 v{<i>) - - Qr x m^i> 

where <& is the dilaton field, <p is an arbitrary scalar field 
and tp is an arbitrary fermion. is the gauge covari- 
ant derivative corresponding to gauge fields with field 
strength F^. k 2 = 8ttG and fi = e ~ K<f,/v/I which is 
the conformal factor used to move from the string frame. 
Powers of ^( < E > ) multiplying terms in the action indicate 
the dilaton dependence of the coupling constants and 
masses. 

More concretely, (f> is Higgs field and V(4>) is its poten- 
tial which we assume to be given by hand (as is done in 
the Standard Model). The all over factor before the 
scalar potential means that the Higgs vacuum expecta- 
tion value (H) is independent of the dilaton so has same 
value for now and BBN time. (H) is taken to be constant 
in our calculation. 

F^ v is the gauge field with gauge group including 
SU(3) x SU(2) x U(l). We define its Lagrangian density 
as —-^F ilv F^ u where g is the unified coupling constant. 
Comparing with Eq. (|l5|), 

1 _ a'e"^* _ a'S 
g 2 (M 2 ) ~ 4^2 (1 ' 

where we define 

S=- . (17) 

K 

For each value of S, we can calculate the gauge coupling 
constants at the low energy using renormalization group 
equations, cxem almost does not run so, 
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a BM (M p ) = ^^ = ^ 7 S- 1 . (18) 
4-7T not' 



~>bbn 



Sue 



1 



OLEM,now 

where we define fractional S variation, 

Sbbn Snow 



1 + V S 



v s 



(19) 



(20) 



From the solution of 1-loop RGE for the SU(3) coupling- 
constant 53, of which integration constant is determined 
by S3 (Aqcd) = oo, 



J^_g-i _ g 3 (M p ) 2 



Yin 



An 27\og{M 2 /A 2 QCD ) 



Aqcd = M p exp I - 



2n 2 a'S 



9ft 



(21) 



(22) 



Therefore, we obtain 

A-QCD,bbn 



A 



QCD,now 



exp 



exp 



2n 2 a'[Sbbn — S now ] 



9k 



8tt 2 



9g(M p ) 



TVS 



(23) 



where we use g(M p ) 2 = 0.1. 

Finally, ip's are the ordinary SM leptons and quarks. 
As we take (H) = const., Yukawa couplings y depend on 
dilaton as oc e K */^/2. I n terms of S. 



Vbbn 



1 



(24) 



V. CONSTRAINTS ON THE VARIATION OF 
THE COUPLING CONSTANTS 

Using the model described in the previous section, we 
can express the coupling dep end ence of the BBN input 
parameters considered in Sec. Ill with the fractional vari- 
ation of the dilaton vs defined in Eqs. ( |l7| ) and (p0|). In 
order to quantify how much variation is consistent with 
the observation, we calculate the abundances for differ- 
ent values of vs in addition to r\ with the standard BBN 
code p3|. Then we calculate x 2 {Vi v s) as > 



(a 



th 



a 



obs\2 



5(af)2 + « fcs ) 2 



(25) 



performed 1000 Monte-Carlo simulations using the values 
of Ref. 0] for the nuclear reaction rates uncertainty and 
Particle Data Group Q for the neutron lifetime 885.7 ± 
0.8 s. For the observational errors for D, 4 He and 7 Li, we 
adopt: 

(D/H) ofcs = (3.0 ± 0.4) x 10~ 5 , (26) 
Y obs = 0.238 ± (0.002) stat ± (0.005) sysU (27) 
log 10 [( 7 Li/H) obs ] = -9.76 ± (0.012) stat ± {0.05) syst 

± (0.3) o «H. (28) 

Eq. @ is taken from Ref. |l| and Eq. @ from Ref. @ 
where the first error is the statistical uncertainty and 
the second error is the systematic one. Eq. (p8|) is from 
Ref. JTtJ with the error we have added for the uncertainty 
in chemical evolution |isj| . 

The results are shown in Fig. |l|. The shape of the 
contour lines are easily understood. The contours drawn 
from the \ 2 °f three elements (Fig. [j] (a)) are just the 
product set of D (Fig. |l| (b)) and 4 He (Fig. (c)) because 
7 Li (Fig. [j] (d)) does not give much constraint. 

D and Li do not vertically constrain much and 4 He 
alone constrain vs- The reason is D and 7 Li abun- 
dances are determined mainly by charged-particle reac- 
tions while 4 He by freeze-out ratio. The former is de- 
termined by Coulomb barrier penetrability or electro- 
magnetic coupling olem and the latter by neutron-proton 
mass difference A m. A m is depend linearly on Aqc d as 



described in Sec. Ill A . aEM and Aqcd are related to 
each other through Eqs. (|l9|) and (|2l]) which came from 
general notion of the gauge coupling unification. This 
AAqcd _. 80tt 2 Aa E M 



relation tells that 



where 



^QCD ,nouj ^ O^E M, now 

Aqem = o.EM,bbn — ctEM,now e ^c. so small increase in 
olem is accompanied by great increase in Aqcd and Am. 
This is the reason why 4 He constrain vs much more than 
D and 7 Li do in our model. 

To make sure, we explain the trend of the He con- 
tour. As vs increases, Aqcd decreases and hence Am 
increases. This makes ^\/reezeout decrease, so 4 He abun- 
dance decreases. Since the 4 He abundance is an increase 
function of 77, the increase in vs relaxes the constraint for 
higher r\. This shows up in the trend that the contour 
goes up in the direction of increasing rj. 

The constraint on the dilaton field variation is obtained 
from Fig. |Tj (a). For 95% confidence level, 



- 1.5 x 10~ 4 < v s < 6.0 x 10" 



(29) 



or using S = er^ 2 ^ / n with « = {SttG) 1 / 2 = (2.43 
10 18 GeV)-\ we obtain 



1.0 x 10 14 GeV < A$ < 2.6 x 10 14 GeV 



14 



(30) 



where i is the type of the element with which we try to 
put a constraint. To estimate theoretical errors, we have 
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(a) D+ 4 He+ 7 l_i 



(b)D 




FIG. 1. The contour lines on the r\ — vs plane drawn (a) 
using the data of the three light elements D, 4 He and 7 Li, (b) 
using D, (c) using 4 He and (d) using 7 Li. The vertical axis is 
Vs _ s bb „-s noqll ^ rp ng gouc ] U ne j g connc i ence level and 

the dotted"line 70%. 



D+ 4 He+ 7 Li 



s » 




1U 1U 

FIG. 2. The contour lines when only ccem is varied drawn 
using the data of the three light elements D, 4 He and 7 Li . 
The relation v s » - AaEM holds. 



VI. CONCLUSION 

In summary, we have considered the BBN calculation 
with the coupling constants taking different values from 
the present ones and whose variation is determined by 
the dilaton field <E> as in the action (|l5|). The set up is 
the same as Ref. Q but we have performed full nuclear 
reaction network calculation with the BBN code so the 



abundances of D and 7 Li in addition to 4 He are avail- 
able to compare with the observation and make a con- 
straint. This is done by calculating \ 2 as a function of the 
baryon-to-photon ratio ij and the fractional variation of 
the dilaton vs defined in Ens. ( p~7[ ) and (p0|), and we have 
obtained the constraint (|29|), which is the same order as 
the value found in Ref. M7 It might be interesting that 
the standard BBN («s=0) is outside the allowed region 
with 70% confidence level as seen in Fig. [l] (a) although 
it is not statistically significant. 

The constraint on vs is readily converted to that on 
olem using Eq.([l^), which becomes Aolem /olem = —vs 
when vs is small, and gives 



_ 6 .0xl0- 4 <^^<1.5xl0- 4 , 
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(31) 



which is two orders stricter than what found in Ref. j^] 
where only olem is varied. This feature should be com- 
mon in the models in which variation of the gauge cou- 
plings are related to each other through gauge unifica- 
tion so the variation in Aqcd is much larger than that 
of olem ■ 

For comparison, we draw a contour diagram when 
only olem is varied. The result is Fig. ||. The shape 
is similar but noticing the vertical scale, it is two or- 
ders larger than the Fig. [j] (a). The constraint becomes 
-5.0 x 10~ 2 < &a em, only I 'ol E m, only < 1.0 x 10~ 2 and 
this is similar to the result obtained in Ref. j5j (factor 
difference is attributed to the difference in the adopted 
observational data). The contours in the Fig. ^ com- 
pletely covers the contours in the Fig. [I] (a). This is an- 
other demonstration that the constraint is not obtained 
by the variation of olem directly but by Kqcd- 

Finally, we notice that the limit (|3l]) is consistent 
and has favourable trend with the result found in 
Ref. § from the quasar absorption systems, that is, 

^OLEM.onlyloLEM.only = (~0.72 ± 0.18) X 10~ 5 . This 

shows olem was smaller at z = O(l). So if the variation 
is monotone, it should be further smaller at the BBN 
time. This is just the tendency of Eq. (|3l]) (the allowed 
region is wider for the negative side). More quantitative 
comparison needs a model for time evolution of the cou- 
pling constants (for the model considered here, it is spec- 
ifying dilaton potential) and the analysis of the quasar 
data under the situation where the couplings other than 
olem vary. 

We have not reached a stage where such a constraint as 
( |30| ) can be used to say something about particle physics. 
In this paper, we bear string theory in mind and treat the 
time variation of the coupling constants as caused by the 
dilaton, scalar field found in the theory. However, in cos- 
mology, there flourishes scalar fields to attack important 
problems (e.g. inflation, quintessence) so the similar anal- 
ysis (with different actions) can be used to constrain such 
scenarios. In order to find useful application, the analysis 
itself has to be sharpened because there still remains un- 
satisfactory estimates concerning the nuclear force. This 



G 



sector will be improved with the field theoretic (QCD 
based) understanding of the nuclear force. 
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